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Abstract 

Oh! We study the classical decay of unstable scalar solitons in noncommutative field theory 

^ ■ in 2 + 1 dimensions. This can, but does not have to, be viewed as a toy model for the 

decay of D-branes in string theory. In the limit that the noncommutativity parameter 9 
^ \ is infinite, the gradient term is absent, there are no propagating modes and the soliton 

does not decay at all. If 9 is large, but finite, the rotationally symmetric decay channel 
can be described as a highly excited nonlinear oscillator weakly coupled to a continuum 
of linear modes. This system is closely akin to those studied in the context of discrete 
breathers. We here diagonalize the linear problem and compute the decay rate to first 
order using a version of Fermi's Golden Rule, leaving a more rigorous treatment for 
future work. 

* Dedicated to Rudolf Haag, on the occasion of his eightieth birthday, with admiration 
and affection 
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1 Introduction 

Solitons are interesting objects. Originally discovered in experiments on shallow water 
waves one and a half centuries ago, they have since then become a rich and fascinating 
area of theoretical and mathematical study. Today, solitons are playing an increasingly 
important role in modern field and string theories. 

A general issue in the analysis of solitons is stability, in both the dynamical (small 
perturbations of the initial conditions within a theory), and the structural (small de- 
formations of the theory) sense. Stability can often be established on the basis of 
topological properties of the solution, or by the integrability of the underlying theory. 
On the other hand, a problem that has attracted much less attention concerns the 
dynamical long-time evolution of decaying, unstable solitons. The reason is apparently 
not because this question fails to be interesting, but rather because it is hard. 

In this paper, we study the decay dynamics of a particularly simple class of sohtons 
that were discovered only very recently, coined noncommutative solitons [201 • The 
corresponding model can be derived from the action functional 



which describes a real scalar field ^{x,y,t) in 2 + 1 dimensions. Noncommutativity is 
introduced in the multiplication rule for $, as we review below. We will only consider 
space-space noncommutativity in this paper, governed by a single noncommutativity 
parameter 6 G [0, oo). 

"Noncommutative field theories" of this kind appear for instance as world-volume 
theories of D-branes in string theory by turning on a large Neveu-Schwarz-Neveu- 
Schwarz B- field. In general, the resulting theory is more complicated than (Q), and 
will include gauge fields and so on, but here, we will restrict ourselves to (^. We will 
not indulge here in a general discussion of noncommutative geometry IJ, and its role 
in modern physics. As starting points for references and citations we refer to the short 
hst PElillSlini, as well as to the reviews U\EB. 

1.1 Motivation from string theory 

Our interest in this problem arose from recent studies of tachyon condensation in string 
theory, hence let us briefly review this background motivation here. Tachyon condensa- 
tion in string theory is formalized in Sen's conjectures [TlU lllj. Roughly, given a string 
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theory background with an unstable D-brane or D-brane system, i.e., given a solution 
of string theory with an open string tachyon, one should look for another solution of 
string theory that is related to the first solution by "tachyon condensation" . This solu- 
tion corresponds to a different background with other branes, and can either be stable 
or contain other tachyons. In bosonic string theory, all Dp-branes are unstable, and 
open string theory is conjectured to have a solution corresponding to the condensation 
of all these tachyons, characterized by: 

(51) There are no perturbative open string excitations around this solution, which is 
hence also known as the "closed string vacuum of open string field theory" . 

(52) The energy difference between the solution with open strings and the closed string 
vacuum is exactly equal to the D-brane tension computed in perturbative string 
theory. 

(To avoid confusion, we note that the bosonic string in 26 dimensions also has a tachyon 
in the closed string sector, whose role is, however, much less understood.) 

A useful analogy in this context is to think of D-branes as "tachyonic solitons" . For 
example, a codimension-one brane in bosonic string theory is analogous to the solitonic 
lump solution in a scalar (p^ theory. This soliton and the D-brane are of course unstable. 
In superstring theory, a codimension-one brane is analogous to the stable solitonic kink 
in a scalar field theory with a (p^ double-well potential. The topological charge of such a 
soliton that guarantees stability is analogous to the Ramond-Ramond charge of a brane. 
These analogies break down at some level because string field theory contains many 
more fields than just the tachyon, and they are all excited in the vacua representing 
the endpoint of tachyon condensation. 

Tachyon condensation can be formulated on the worldsheet (as worldsheet RG flow) 
or in spacetime in open string field theory ^2]- This has been quite successful and can 
be regarded as one of the major advances in the field over the last few years. Until 
recently, however, the research has focused on finding the endpoint of the condensation, 
and not on the decay of the unstable brane (the condensation of the tachyon) as a 
dynamical process in time. 

Very recently. Sen p!H | ll4 | [TK] has made precise proposals for the construction of 
the time-dependent solutions corresponding to the condensation process. This is based 
mainly on a worldsheet analysis, and is motivated in part by the S-brane proposal of 
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Gutperle and Strominger ^Hl- Again, there are many different possible processes one 
can study. The generic features characterizing the solutions are 

(53) If the initial configuration is an unstable brane configuration of non-zero codi- 
mension (say a Dp-brane in bosonic string theory with p < 25) then the energy 
stays localized in the plane of the original brane, in other words, energy is not 
carried away from this plane by classical radiation. 

(54) When computing the energy-momentum tensor of the solution on finds that, 
asymptotically for late times, the energy density of the solution is constant (in 
the plane of the original brane) and the pressure vanishes exponentially in time. 

It is important to emphasize that these are supposed to be classical solutions of 
string theory. In particular, the radiation in (S3) are open strings. It is very likely 
that in fact closed strings (which appear at one loop level in open string theory) do 
carry away energy [TH t lTT t lTS] . Furthermore, one may note that (S3) is in some sense a 
consequence of (SI). Namely, if one thinks of a lower-dimensional brane as a tachyonic 
soliton, then far away this solution looks like the closed string vacuum which does not 
allow open string excitations. Hence, there are no modes that could carry away energy. 
Also, (S2) and (S4) are of course partly trivial consequences of energy conservation. 

All these questions are rather difficult to study in open string field theory. The 
underlying idea behind the present work was to study the dynamics of such decays in 
the so-called noncommutative (NC) field theory limit of string theory, in the hope that 
one could reproduce some of the above properties in a simple manner. 

The reason for this hope is the fact that solitons in noncommutative field theories 
(NC solitons) have many properties in common with D-branes as tachyonic solitons. 
NC solitons and instantons in gauge theories have been studied independently of strings 
for quite some time, most notably by Nekrasov and Schwarz JHI; and more recently 
by Harvey et al. [21], Gross and Nekrasov |221I211; Polychronakos and others. 
In NC scalar field theories, solitons were discovered by Gopakumar, Minwalla and 
Strominger [20] • For a summary of the relation to D-branes, in particular as far as 
their charges and descent relations are concerned, see Harvey's lectures [2E\ . 

An explanation of why NC solitons (in the infinite 6 limit) are useful for describing 
tachyon condensation in string theory was given by Witten [2B]. The point of j2S] 
is that in the Seiberg- Witten limit, the open string field algebra factorizes into the 
Moyal algebra built on tachyon vertex operators and another algebra that "contains 
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all the stringy mystifications" . Therefore, in describing tachyon condensation, one can 
concentrate on the tachyon, the remaining pieces being added "trivially" at the end. 

Given all this, it becomes a natural problem, from the point of view of string theory, 
to study the decay of unstable NC solitons, using (S1)-(S4) as guiding questions, as a 
toy model for the dynamics of open string tachyons on unstable D-branes. 

1.2 Summary of Results 

We will see that it is in fact not possible to reproduce the qualitative behavior described 
above in the simple model ((H). A priori, one would not expect (S3) simply because 
of the failure of (SI) in the context of NC tachyon condensation. Indeed, while NC 
solitons have localized energy, far away the field is in the usual vacuum that has local 
propagating excitations. However, we note that intuitively (SI) is sufficient but not 
necessary for (S3), and in this sense, finding (S3) would have been a nice substitute for 
(SI). And recall that the very existence of NC solitons at finite 9 relies on the fact that 
the kinetic term is a small perturbation that does not destroy the qualitative properties 
at infinite 9. One could expect something similar in the time-dependent case. 

What we will find is that, except for infinite 9, the energy that initially is stored 
in the soliton does not stay localized, but is carried away to infinity in the form of 
classical radiation. Nevertheless, the qualitative properties of the decay do depend 
on the noncommutativity parameter 9. In particular, the decay rate is exponentially 
suppressed in 6*, for large 9. 

We will proceed as follows. We start with a review of NC scalar solitons and their 
basic properties in section El It turns out that restricted to rotationally symmetric 
configurations, the system is similar to certain nonlinear lattices which exhibit the 
phenomenon of "discrete breathers", and we review those in section 01 for completeness. 
We will then present some numerical results in section ^ in order to illustrate the 
qualitative properties of the decay. In section we turn to reproducing this behavior 
through an approximate calculation of the decay rate. The technical core is here the 
explicit diagonalization of the linear problem in section (HI that puts us in a position to 
give a closed expression for the estimated decay rate in section [3 The approximation 
method relies on a certain version of Fermi's Golden Rule that we will explain. 

To conclude this introduction, we note that a different toy model for tachyon con- 
densation in string theory, based on conventional field theory, has been considered by 
Zwiebach |2H1- It would be interesting to compare the dynamics of the decay with our 
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present results. Also, Moeller and Zwiebach have studied the time evolution of 
decaying solitons in p-adic string theory. The relation to our work is not clear. Finally, 
we note that some of the problems that we study here have also been broached in [HI] 
and |H2|. The stability of non-commutative scalar solitons is also discussed in p?!^l34j . 



2 Noncommutative scalar field theory in 2 + 1 dimensions and 
its solitons 

The theory of our interest is defined over a model of 2 + 1-dimensional spacetime 
in which time is an ordinary variable and space is the two-dimensional quantum me- 
chanical phase space of a nonrelativistic particle in one dimension. In other words, we 
base space on the Heisenberg-Born- Jordan commutation relation, 

[x,y] = te. (2) 

It is well-known since the times of Weyl, Wigner, and Moyal that the algebra 
of functions on this space is a deformation of the ordinary algebra C(M^) along the 
parameter 9. In one particular ordering prescription, the product can be written as 

$1 = e'i(^^i^--^^^^"i)<l>i(6,r/i)<f2(6,^2)|(6,.i)=(6,..)=(x,.) • (3) 

Upon representing the commutation relations Q in the usual manner on Hilbert 
space 7i, eq. turns into a prescription for translating a function over phase space, 
such as our field $, into a linear operator on Ti. Moreover, derivatives become com- 
mutators, 

dx< — >Jy,-], dy< — '-nl^,-], (4) 



and the integral a trace. 
This yields, after setting 



dxdy< — i>27r^tr. (5) 



X + iy) , = [x — iy) , (6) 



^29 V29 
the following expression for S. 

S = 27i9 j f/ttr[^(9i$)2-i[a,$][$,at]-y($)] (7) 

We also record the equations of motion 

a,^$ + ^[a,[at,$]] + V^'($)=0. (8) 
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2.1 Rotational symmetry 

These expressions simplify further on rotationally symmetric field configurations, 

oo 

$(t) = ^A„,(t)P„. (9) 

n=0 

Here, Pn = \n) (n| is the projector onto the n-th excited state in the harmonic oscillator 
basis for H. One finds 

[a, [a\ Pn]] = {2n + 1)P„ - (n + l)P„+i - nP^-i . (10) 

The equations of motion (jHJ become 

2 

Ao = ^(Ai-Ao)-y'(Ao). 



'11^ 



This is a hamiltonian system with Hamiltonian 

oo ^ 

= 27r^^ 5^ [-kI + ^(A. - A„_i)' + ^(An)] , (12) 

n=0 

where the momenta canonically conjugate to the A„. 

To connect (jllj) and ()12|) to something more familiar, we note that using Q to 
write the projector P„ as a function on M?, it has a peak around 

r'^ = {n\x^ + y^\n) K. 2dn , (13) 

and that substituting ()13|) into the radial part of the Laplace operator written in two- 
dimensional polar coordinates, one obtains 

-dr-rdr = \dnndn . (14) 
Upon discretizing n and writing dn as a difference operator, one reproduces (|TT|) or 

m- 

In this paper, we study properties of the dynamical system (fTTj). 
2.2 Non-commutative solitons 

We start with static solutions of (jH)) and ()11|) . the noncommutative solitons. It was 
pointed out in [20] that in the limit 6* — cx), eq. (jHl) simply reduces to V'{^) = and 
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can be solved by a projector $ = XP with A G M such that V'{X) = 0. One then 
treats the kinetic term as a perturbation and shows that the solutions continue to exist 
for 6 < oo. Since, according to Derrick's theorem, the solutions do not exist in the 
classical or "continuum" limit 6—^0, there must be a critical 6c beyond which the NC 
solitons cease to exist. Of course, existence is not controlled by 6 alone, but rather by 
certain dimensionless quantities built out of 6 and the parameters of V. 

The NC scalar solitons which are related to projectors in the 6 oo limit are re- 
ferred to as the GMS (Gopakumar-Minwalla-Strominger) solitons. Instead of reviewing 
the large number of works devoted to GMS solitons, let us quote an existence theorem 
proved in [37], see also pUj . 

Theorem 1. For every sequence X^^^ = {X^q\ Xf'\ . . .} with V'{Xn^) = and 
V"{Xn^) > for all n (i.e., a stable static solution of (fTT|) at 1/6 = 0), there ex- 
ists a unique continuation X{6) = {Ao(^), Ai(^), . . .} to static solutions of (fTTj) for 

< 1/6 < l/6c, with some 6^ which depends on the chosen initial sequence 

We sketch the proof given in YJlj . We work on the Hilbert space /2(^>o) = Ti. = 
{^'^Y^n=o\^n\'^ < Let {AX)n = -{u + l)(A„+i + A„) - n{Xn - Xn^i) be the 

discretized radial part of the Laplacian according to (jl4|l . with domain T). 
The equation we have to solve is then 

eAX + V\X)=Q, (15) 

where e = 2/6. 

The first step of the proof is to show that the operator eA + V"{X'^'^'>) is invertible 
with bounded inverse as a map from T> to Ti. This essentially follows from l^"(A*^°^) > 0. 
One then rewrites as 

X={tA + V"(A(°)))"' [r"(A(°))A(°) + \/'(A(°)) - V'{X) - V"(A(°))(A(°) - A)] =: r,(A) , 

(16) 

and shows that : ^ 7i is a contraction in a neighborhood of A*^°^ (for reason- 
able regularity assumptions on V). The Banach fixed point theorem then implies the 
existence of a solution, with smooth dependence on e in the appropriate norms. □ 

2.3 Stable and unstable solitons 

In order to simplify the following discussion, we assume that the potential V is such 
that it has one global minimum Xmin, one local minimum Amm 7^ Xmin, and one local 
maximum Xmax- Corresponding statements should hold for a general potential. 
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In ref. [37], certain statements are proved concerning the stability of the solutions 
whose existence for 1/9 > is established by the above Theorem ^ Stability is here 
defined as positivity of the Hessian of the Hamiltonian = 0, A), eq. (jE}, viewed 

as a quadratic form. For example, if V^(A) is as decribed above, then the solutions 
obtained as continuations of A^^-* to 1/^^ > are stable if and only if A'-"-* is of the form 
Aq°^ = A^"'* = ■ ■ ■ = Aj^^ = Xmin for some N and Xn^ = Xmin for n > N. 

It is interesting to note that the unstable direction of the unstable solitons of The- 
orem [T] [e.g., the one connected to Aq°^ = Xf^ = ■■■ = Xmin, Xf'^ = Xmin) is not 
rotationally symmetric. More precisely, there exists a non-rotationally symmetric in- 
stability whenever the sequence Ao(6'), Ai(^), . . . at finite 6 is not monotonic. The 
nature of the instability implies that the decay of these solitons can not be studied 
with an ansatz of the form Q. 

Here, we shall be interested in a different kind of unstable soliton, namely the one 
obtained from continuation to 1/^^ > of the static unstable solution Aq'''' = X^ax, 
Xn'^ = Xmin for n > 0. It is in fact rather simple to adapt the above proof to show 
existence of these solitons as static solutions.^ The essential point is that eA + V"{X^^^) 
is still invertible even if r'(AS')) is negative for a finite number of n's, but positive 
elsewhere. This is because it differs from a positive operator only by a finite rank 
perturbation, so that eA + V"{X^'^'^) can still be inverted for e in a finite neighborhood of 
0. This soliton is unstable, but the instability is rotationally symmetric. Moreover, the 
sequence of A„'s is monotonic, so that at least initially, the soliton does not suffer from 
the non-rotationally symmetric instability discussed above. This justifies an ansatz of 
the form (jU)). 

Summarizing, we can imagine studying the decay and long time evolution of the 
following unstable initial conditions, slightly perturbing them away from the stationary 
point. 

(11) A GMS soliton at finite 6 obtained from an unstable static solution at infinite 6, 
say with Aq"^ = Xmax, An ^ = Xmin for n > 0. 

(12) A GMS soliton at finite 6 obtained from a stable configuration at infinite 6, say 

, (0) _ , (0) _ _ , ,{0) _ , 

^0 — ^2 — • ■ ■ — AMIN, ^1 — ^min- 

(13) The unstable maximum, A„ = Xmax for all n. 
^We thank M. Salmhofer for this suggestfon. 
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The questions to ask about the decay are quite different from case to case. In view 
of the relation between the noncommutative tachyon and string field theory explained 
in j2ni, one can draw an analogy between (II) and the decay of an unstable D23 brane 
in bosonic string theory, in which one looks at the directions transverse to the brane. 
In view of (S3), then, one can ask if and how fast the energy that initially is localized 
is transported away to infinity. We will see that while energy is indeed dissipating 
whenever 1/6 > 0, the decay is very slow. 

We emphasize, again, that the action is not a good model for the tachyon of 
string field theory, even in the large B-field limit. Rather, as explained in [25], not only 
does the algebra factorize, but also the equations of motion do, and the tachyon has 
no dynamics in this picture. Thus, in the large B-field limit, one always obtains the 
noncommutative "field" theory at infinite 6. It is a nice, but rather trivial observation 
of consistency that in this limit the energy does stay localized j31j . 

We find it harder to give an interpretation of the decays that start from (12). We will 
not study these decays here, but we expect that some generalization of the technology 
that we develop could be applied in that case. 

We finally comment on the initial conditions (13). Again by analogy with string 
field theory, one might regard (13) simply as the decay of the open string vacuum (or 
of a lower dimensional D-brane in which one looks in the directions along the brane). 
Then, one could ask for comparison with property (S4) of open string tachyons, in other 
words, can one see the formation of "tachyon matter" 11211111113 • We will not study 
this here, by let us observe that if the perturbation away from the stationary point is 
uniform, i.e., independent of n, then the equations of motion can be trivially solved to 
yield uniform oscillations. So, starting from (13), one may want to ask what happens 
under non-uniform perturbations. It also makes little sense to restrict to rotationally 
symmetric decays in this case. 

In this paper, we will only deal with initial conditions of the form (II). 

3 Discrete breathers 

In order to explain some intuition about the behavior of the system we here 

digress on the theory of discrete breathers. 

Discrete breathers are localized periodic solutions of lattices of oscillators, in which 
localization is not due to disorder (as in the usual case of Anderson localization). 
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but to nonlinearities. A useful review is |]^. The original discovery can be traced 
back to [SHlEni! in which numerical simulations on these lattices were performed. A 
rigorous existence proof was given in PUI, further and related earlier results appear 

in nniiaiisiiiaiisi. 

As an example, let us consider a Hamiltonian similar to (fT2|) . 

oo ^ 

H= + V{Xn) + e„(A„, - Xn-lY] . (17) 

?i=— oo 

For en = 0, this is just an infinite collection of independent classical oscillators. For 
oscillations around a minimum A min of V (we have in mind a potential similar to the 
one described in the previous section), we introduce canonical action-angle variables 
(J„,<^„), with 

Jn = ^f V2(^( J„) - F(A„)) rfA„ , (18) 

and we let 0„ = correspond to the turning point of the oscillation, k„ = 0. These 
coordinates are valid until uj{Jn) = ^^j^ = 0. Such oscillations are localized at 
whichever site of the lattice is excited. 

On the other hand, if we couple the oscillators with e„ = e 7^ but small and 
constant, the lattice has propagating small oscillations (the phonons). The spectrum 
of the linearized problem is 

u;lf^{k) = + 4e(l - cos k) , (19) 

where uJq = V"{Xmin), and k G [0,27r]. 

The linear spectrum is therefore bounded and given a nonlinear oscillation of the 
problem at e = that is not in resonance with ujq, i-e., muj{Jn) 7^ ojq for all m & one 
can avoid resonances with the phonons, muj{Jn) 7^ ujph{k) for all m G Z and k G [0, 27r], 
if e is sufficiently small. It is therefore plausible that such an oscillation can continue to 
exist and be localized for e 7^ 0. This intuitive argument can be turned into a rigorous 
existence theorem for discrete breathers [10] . 

Theorem 2. Every periodic solution of the decoupled nonlinear problem with frequency 
uj satisfying the nonresonance condition muj 7^ ujq for allm E 7^ and the anharmonicity 
condition ^^i^^i^l ^ for all n G Z and with initial conditions coded by a sequence 
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Vn e {0,1}, 




ifVn = -1 
if Vn = 0, 



(20) 



can be uniquely continued to a periodic solution of the full nonlinear problem e„ = e, for 
e > sufficiently small. These oscillations are localized in the sense that |A„(t) — Xmin\ 
decays exponentially in n whenever ?7„ — >• 0. 

The proof, which can be found in uses an imphcit function theorem and is 
rather elementary. A more geometric formulation is given in The periodic solu- 
tions given by the theorem are dynamically stable, in the sense that small perturbations 
of the initial conditions remain small for all times. This is discussed in an abstract 
setting in . 

We note that the solution that one wishes to continue has to be strictly periodic at 
e = 0, i.e., all excited oscillators have to be in resonance and they also have to be in 
phase, as coded by pO|l . This resonance condition appears to be opposite to the one 
found in the formulation of the KAM theorem. However, it is completely in agreement 
with the intuitive argument. If two frequencies are not rationally related, there will 
always be a resonance of some higher harmonics with the phonon band, leading to 
decay. 

What can discrete breathers teach us about NC soliton decay? 

The analogy with discrete breathers leads to the question whether breather-like solu- 
tions could exist in NC field theories. The answer seems to be no. 

The obvious difference is that in ()12|) the coupling en = n/6 is not bounded as soon 
as 6* < oo. The linear spectrum, which will be discussed in detail in the next section, is 
therefore unbounded. This leads to a failure of the intuitive argument given above for 
existence of breathers. Any nonlinear oscillation will be in resonance with the phonons, 
no matter how small 1/6. There could be some hope if the linear spectrum was pure 
point, as in Anderson localization. However, as we will see below, the spectrum is 
purely absolutely continuous. 

It is a little more difficult to understand what goes wrong with the proof of |40j . 
But recall that one of the assumptions for the implicit function theorem is continuity of 
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the derivative, and this seems to fail here, because the derivative with respect to 1/^ is 
unbounded. While there might still be some hope that one could use slightly stronger 
versions of the implicit function theorem to prove existence, this will generically be 
restricted to very special conditions, and one cannot expect these breathers to be 
stable. In any case, our numerical results and the approximation computation of the 
decay rate clearly disfavor the existence of "noncommutative breathers" . 



4 Numerical Results 

We now turn to presenting some results obtained by studying the system (jllll using 
numerics on the computer. We use the potential 

=4$2($-2)2 (21) 

and initial conditions of the form (II). Specifically, we take 

-^nlO) = Amor ~ 1 

^ ' A„(0) = (for all n) . (22) 

A„(0) = Xmin = (n > 0) 

For nonzero 1/6, this is not a stationary point, and will lead to some decay. 

We have to work with a finite system < n < A^. Therefore, in order to see a decay, 
if it occurs, we have to allow for a leak of energy. We do this by imposing "absorbing 
boundary conditions", i.e., the discrete version of 

^ + ^ ^ Vr-0 (23) 
dt yfr dr ' 

(recall that outgoing spherical waves in two dimensions behave asymptotically like 
g«E(t-r)^^y^^^ which in view of (fT!^ reads 



A„« + ^?^^(A.«-A.) + ^=L=A„,.=0. (24) 

We start by plotting, as a function of time, the energy Ei^c that is stored in the 
first four (say) oscillators, i.e., 

Eiocit) = Y.[\\l + ViK) + ^(A„ - A„„i)^] , (25) 

71=0 

divided by the initial energy, Eioc{Q). For different values of 6, our results are shown 
in figs, [mi 
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Figure 1: Energy ()25j) stored in the first four oscillators of the system flllj] . truncated 
to 21 oscillators with boundary conditions (A^ = 20). The initial conditions are 
dUl) forn < and Aiv+i(0) = 0, and 9 = 0.1. 
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Figure 3: Same as fig. Q for = 0.35. 
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Figure 4: Same as fig. HJ for 6^ = 0.7. 
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While the energy seems to always dissipate, there are apparently significant differ- 
ences between large and small 6. For small 6', the system dissipates most of the energy 
at early times. For large 6', the system loses a small fraction of its initial energy rather 
rapidly at early times, and after a certain "transient" time the energy dissipation rate 
approaches a constant but rather small value. To see this more clearly, we estimate 
the energy that is left after the transient time by averaging Eioc{t) over a certain late 
time interval, and plot the result as a function of 6. The result is shown in fig. El 




Figure 5: The energy left in the first four oscillators after the initial transient time 
apparent in figs. miH obtained by averaging Eioc{t)/ Eioc{0) from t = 30 to t = 50, for 
varying 6. 



We see that the turnover happens roughly at 9 ^ 0.23. This agrees rather well with 
the estimated critical 6c below which the unstable NC soliton is expected to disappear. 
From the most naive estimate is 



0. ' 



V'iX)\ 



^ 0.26 (26) 

V"(\)=0 



for the potential (j?T|) . We also note for comparison that the period of linear oscillations 
is To = 27r/cuo ^ 1.1, where = V"{0). 

In the remainder of the paper, we will compute the rate of energy dissipation to a 
first approximation in the coupling 1/6. 
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5 The Setup 



We start from the classical Hamilton function 

oo _ 

H = Y. \r/n + ^(A.) + ^(An - K-.f\ , (27) 

n=0 

describing an infinite collection of oscillators, labeled by n = 0, 1, . . ., with n-dependent 
coupling. We want to understand the late time behavior at large ^, if at t = we 
excite the zeroth oscillator highly into the non-linear regime. The numerics indicate 
that all other oscillators remain weakly excited during the decay process, so that we 
may linearize around the minimum of V . Furthermore, for large ^, the zeroth oscillator 
is weakly coupled to all the other ones (the other oscillators are strongly coupled as 
n oo). Accordingly, we split the Hamilton into three pieces 

H = Hq + Hcont + -f^coup , (28) 

where Hq describes the zeroth oscillator that we want to excite, 

^0 = \4 + ^(Ao) , (29) 

-f^cont describes the remaining oscillators linearized around the minimum, which for 
convenience we assume to be at \min = 0, 

Er 2 2 \ 2 ^ I / \ \ \ 2~l / \ 

V^i^n + 2^oA„ H — - A„) J , (30) 

n=l 

and c<Jq = V"\min- Finally, ifcoup describes the coupling between Hq and -ffcont 

-f^coup = ^(Ai — Ao)^ . (31) 
At this level, the equations of motion are 

Ao + ^'(Ao) = ^(Ai-Ao) 
Ai + ^o'Ai - ^ [2(A2 - Ai)] = -^(Ai - Ao) (32) 

2 

An + u;o A„ - - [(n + 1) (A„+i - A„) - r2(A„ - A„_i)] = forn > 2 

Our strategy is as follows. We intend to treat the RHS of the first two equations in 
as a perturbation. Without this term, the first line can be trivially integrated. By 
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our assumptions on the potential V, the solution, which we call Aq, will be periodic, 
and we call the corresponding frequency a = where Jq is the action variable. 

We then feed this oscillation into the last two equations of (jS2I), which we write as 

2 2 ~ 

A„ + cUgAn, + -(AA)„ = -5„,iAo , (33) 

where A is the linear operator 

{AX)n = -{n + l)(A„+i - A„) + n(A„ - A„_i) , (34) 

supplemented with appropriate boundary conditions (see below). We can view (jHHj) 
as a radiation problem with Aq as the source. It is then natural to decompose Aq into 
Fourier modes, 

oo 

Ao(t)= (35) 

m=— oo 

and to first compute the energy radiated into the continuum modes of A by a harmonic 
source. Summing over all Fourier modes, this will give an estimate for the rate of 
energy loss of the zeroth oscillator due to the coupling to the continuum modes. The 
expectation is that for large 6, the energy loss per period is small, so that the coupling 
to the continuum essentially amounts to adding a friction term for Hq , with dissipation 
rate depending on the remaining energy. One can in principle study the resulting one- 
dimensional dissipative system, but we have not done this here. On the other hand, 
for small 6, in particular 6 < 6c, the energy loss per period will become comparable to 
Ho . In fact, this will give an estimate for the critical 6c- 

Most of what follows is concerned with the operator A. In the following section 
El we will diagonalize A, which means finding the generalized eigenfunctions and in 
particular their asymptotics aX n —>■ oo. This will allow in section [7| to compute the 
power radiated into the continuum by a harmonic source. Finally, we estimate the 
Fourier coefficients of Aq to obtain a decay rate of the soliton. 

6 The operator A 

The operator A defined by (j34p has the form of a discrete version of the familiar class 
of Sturm-Liouville operators, and we will draw on some of the theory pertaining to the 
latter (see, for instance, jlHl)- First, however, we will obtain some elementary results 
on the nature of the spectrum of A. 
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6.1 The spectrum of A is unbounded and absolutely continuous 

We have 

Proposition 1. The spectrum of A is equal to [0, oo) and purely absolutely continuous. 

To prove the first part of the statement, we recall that the spectrum of Aj^, 

(AjvA), = iV(2A„ - A„+i - A„_i) (36) 

on hC^) is equal to [0, 4A^]. The spectrum is absolutely continuous and the generalized 
eigenvectors are A„ = e*'^" with spectral parameter k = 2N{1 — cos A;). 

Now, for fixed k G [0, oo), we consider G Z>o large enough such that k = 
2N{1 — cos A;) has a real solution k{N) > 0. For given N,M G Z, we can find a 
normalized 'almost-eigenvector' A = (A„) G /2(Z>o) with support in [A^, + M] and 

||(A;v - n)X\\ < const. + ^) , (37) 

where the constant is independent of and M. (For instance, we can take A„ = 
f{n)e^''^^^'^, where / is normalized, has support in [A^, A^ + M], and is sufficiently 
smooth.) We also have an estimate on the first derivative A', (A% = A„ — A„+i, e.g., 

\\\'\\ < const. (^k{N) + (38) 

Then, using 

n — N 

{{A - An)X)^ = ^^(A^A)„ + A„ - A^„+i , (39) 

we find 

IMk ,,,,, 1 /M \/Nk(N) N \i 
IK^ - MAII < const [— + k(N) + -g + (-^ + l) (^ + ■ (40) 

For large A^, we have k{N) ~ (k/A^)^/^ and we let M ~ N^/'^. Thus, 

||(A - k)\\\ < const. [kN-'/^ + k'I^N-'I'' + Ar-3/4 ^ ^1/2^-1/4 ^ ^-1/2] ^^^^ 

Making A^ large, the Weyl criterion implies that n is in the spectrum of A. Positivity 
of A is trivial. 

To prove absolute continuity, we use the so-called positive commutator method 
or Mourre theory [47j. The point of Mourre theory is that one can prove absolute 
continuity of parts of the spectrum of a self-adjoint operator such as A by simply 
exhibiting another self-adjoint operator K such that the commutator i[A, K] is positive 
on the part of the spectrum in question. 
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A positive commutator. We write 



A = 2n + l-nT+- T_n , (42) 

where (nA)„ = nA„, (T_|_A)„ = A„_i, and (T_A)„ = An+i. We have 

[n,T±]=±T±, [T+,T_] = -5„,o. (43) 

Define K = |(nT+ — T_n) and find 

ir] = -[2n + 1 - nT+ - T_n, i(nT+ - T_n)] 
= — n[n, T_|_] + [r?,,T_]?7, 

- l{nT4T+,n] + [n,T_]nT+ - [T_,n]T+n - T^n[n,T+]) (44) 
= —nT^ — T^n — | (— 'n,T_T+ — T^nT^ — T^T^n — T^nT^^ 
= 2n + 1 - nT+ - T^n = A > . 

This shows absolute continuity away from 0. But is not an eigenvalue, because any 
solution of n(A„ — Xn+i) — {n + l)(A„+i — A„) = is constant and either trivial or 
nonnormalizable. □ 

We now turn to a more formal analysis of the discrete Sturm-Liouville operator A. 



6.2 Boundary Conditions 

A is naively a second order difference operator, and it might seem that we need to 
specify boundary conditions at = 0. However, from the definition (jSH), it is obvious 
that n = is a singular point of the difference operator (the order of the operator 
jumps there). Therefore, we do not need to specify boundary conditions, if we consider 
the operator on Z2({0, 1, . . .}). We will denote this operator (called A in the previ- 
ous subsection) acting on Z2(^>o) as in (j34|) without boundary conditions, by Aq. An 
alternative definition of Aq would be by a quadratic form. 

For our radiation problem (jH^ . we want to treat Aq as an inhomogeneous 

term added to a homogeneous linear equation (the LHS of (j^Sl))- In order for (jHHj) 
to coincide with (jH^ . we have to define A acting in /2({1,2, . . .}) as in (jSlj), with in 
addition the boundary condition 

Ai - Ao = . (45) 
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This boundary condition makes A self-adjoint. But note that, as we review below, 
the solution of the radiation problem (computation of the resolvent of A) requires the 
knowledge of the solutions of the finite difference equation with arbitrary boundary 
conditions. Also note that the particular form of the boundary condition is a 
result of the way we have separated the perturbation in (jH^ . We could, for instance, 
also have left |Ai on the LHS of the equations. Our results should not depend on this 
choice. 

6.3 The generalized eigenfunctions 

We know already that the spectrum of Aq is [0, oo) and absolutely continuous. The 
same statement holds obviously for A. Here, we want to complement this result by 
explicitly diagonalizing A. From the spectral theory of Sturm-Liouville operators, it is 
well-known that we will need to know the two linearly independent solutions of 

(AA)„ = fi:A„ (46) 

for arbitrary k G C. These two solutions will not necessarily be I2 normalizable and 
satisfy linearly independent boundary conditions at n = 1. Specifically, we call the two 
solutions of (Uni), An = $k('^) and A„ = "^^in), satisfying the boundary conditions 

$«(0) = 1 $,(!) = 1-K (47) 

^.(0)=0 ^.(1) = 1. (48) 

The reason for choosing ()47|) was that $k then also is a solution of Ao$k = but 
this turns out to be unimportant. The reason for choosing ()48|) as boundary conditions 
for is that in this way, the Wronskian of the finite difference equation, which is given 
in the standard Sturm-Liouville fashion, is equal to one, 

W^(<l>,,^,) = (n+l)(<l>,(n)^.(n + l)-^«(n)<l>.(n + l)) ^ ^ 

(49) 

= const. = 1 . 

We note that the equation has been studied in [22] • It seems, however, that 
one of the functions given in [^2] is incorrect, and moreover, the asymptotics of the 
solutions in the regime where we will need them have not been worked out. Quite a 
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bit of work on 



shows that the two solutions are as follows. 



fc=0 

n 



k\ \k 



A:=0 



'-s=0 



n 



(50) 
(51) 



where Hi, are the "harmonic numbers 



)5 2 



k 

Hk = Y.-. 



(52) 



It is trivial that ()50p and ()5H) are entire functions of k G C. 

We see from fl50p . (j3T|) that ^^(n) is a polynomial in k of order ra, while ^{n) is a 
polynomial of order n — 1 (it is easy to see that the coefficient of (— /t)" vanishes). In 
fact, $k(?t.) is nothing but a Laguerre polynomial, which is also a particular confluent 
hypergeometric function. Namely, 



$«,(n) = L„(/t) = iFi(-n, 1; k) . 
A similar, but more complicated statement can be made for "$,^{11), see below. 



(53) 



6.4 Asymptotics 

We wish to determine the asymptotic behavior of ^^(ri) and ^k(^) for k fixed, and 
n ^ \k\ large, for k in the positive half plane C_|_ :={/«€ C|Re(K) > O}. It was 
noticed before that for large n, the operator A, eq. (jSlI), can be approximated by a 
differential operator, A ^ —-^drrdr, where n = r^/2. Let us also recall that the Bessel 
functions of order zero, 

-xV4)^ 



Jo(x) 



fc=0 



ro(x) = -[jo(x)(ln| + 7)-X] 



k=l 



(-xV4y 



Hk 



(54) 
(55) 



where 7 = lim„^oo(-f^n — Inn) is the Euler-Mascheroni constant, are two linearly inde- 
pendent solutions of Bessel's differential equation 



g" + -g' + g = 0. 

X 



(56) 
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We choose the standard branch for the logarithm in Furthermore, it is easy to see 
that both and the series appearing in are entire functions of a; G C Hence, 
it is clear that Jo{x) and Yq{x) are analytic in C+. 

Comparing ()4(ij) with ()5fj|l . we expect an asymptotic expansion of and 
in terms of the Bessel functions Jo{x) and Yo{x), with x = y/AriK. Here, and in all 
similar cases considered below, we choose the branch of the square root for which the 
signs of Im(a;) and lm{K) are equal. 

6.4-i Asymptotics of ^^{n) 

For $K, such an asymptotic expansion is indeed known (see, e.g., PF). It can be 
derived as follows. Because of ^^in) = Ln^n) satisfies the Laguerre differential 
equation 

Kf" + {l-K)f' + nf = 0, (57) 

where differentiation is with respect to n. If we make the usual substitution /(/t) = 
e'^/'^g{K) with n = x"^ /An in eq. ()57|) . we obtain 

where differentiation is with respect to x. 

For |a;| -C n, which is equivalent to n ^ this reduces to (j^Uj) . up to a perturba- 
tion of order C(^). Thus, in this limit, any solution /(k) of is asymptotic to a 
linear combination of Bessel functions Jq{x) and Yq{x), where the deviation from the 
exact solution is an asymptotic series in powers of The particular linear combina- 
tion at hand can be determined from the behavior at x, k — 0. 

From the expressions for $k(^) in eq. and Jq{x) in eq. it is easy to see 
that 

lim $,„(n) = Jo(x), (59) 

n—^QO 

2 

where k„ := for any fixed value of x G C^. Thus, by what has been stated above, 
it follows that for any fixed k G C4. and n ^ |k| large, 

$4n) ~ e'^/Vo(x(n,K)), (60) 

with x(n, k) = \J Atlk. 
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6.4-2 Asymptotics of ^{n) 

The asymptotics of ^kI?^) for large n ^ \k,\ requires more work. It is natural to expect 
a connection to the second solution of Laguerre's differential equation ()57|1 . Surpris- 
ingly, it turns out that this solution is nowhere to be found in the usual references, 
including [1211131101 • The point is that eq. (j^Tj) is a very special case of the confluent 
hypergeometric equation in which both parameters are integer, and the "logarithmic 
solution" given in the references reduces to the ordinary one in this case. One has to 
go back to jSI] to find that a second solution is given by 



fc=i ■ ^ ^ 



k=n+l 

which one can easily verify by direct substitution, and where we choose the standard 
branch of the logarithm. It is easy to see that both series on the RHS of ()6H) define 
entire functions of k G C 

Before we relate Kn{K,) to "^^in), let us determine its asymptotics for large n 3> 
for K, G C+. We use the same logic as for Ln{n) and first consider the limit n — oo, 
K 0, with = Aku fixed. Defining k„ := |^ , we easily see that 

oo / \l. 

El ^ Th) 

k=i (^') (62) 
~ hYq{x) — Jo{x){hin + 27) for n — >• 00 . 

The same arguments that were used for Ln{K,) then imply that for any k G C+, and 
n ^ large, 

Knin) ~ ne^/^Yoixin, k)) - e^^\Jo{x{n, /t))(lnn + 27) , (63) 



where x{n, k) = ^/Ann. 

To relate ^^^(ri) and Kn{K), we rewrite the sums in (|3T|) in terms of s = k — s and 
k = k — s = s. This yields 



5=0 '-k=0 ' 



k 

I II. \ 

(64) 
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Using this, we claim that for k G C+ and n large 

(— k)^ / n 



s=0 



k=0 



<Cmr'^, (65) 



for a constant C that is uniform in n and Indeed, the LHS of ()()5j) is bounded by 

n n+s 



" lit 



s=0 '-k=n-s+l 



s=0 f=n 



2n 



(66) 



<Cr^(2|«:|)"^5:2^ 



since Hi < 2 logn for all 1 < / < n, and (|) < 2". This straightforwardly implies ()65|) . 

By noticing that the last term in (j6ip is small in the limit n ^ \k\, with absolute 
value bounded by 



E 



n! (r — V, 



^ in ^ r\\ in ^ r\\ n\ 

r=l ^ ' ^ ' 

one then infers that for any k G C+, and large n ^ |k| , 



(67) 



E 

•fc=0 



— K / n 



k\ \k 



n—k 



2 Hi, 



E 

s=0 



Si 



-H. 



(— n 



fc=0 



(68) 



~ e-« [K„(/€) + LMiHn - In /t)] + LM ^ ^^^/Z, , 

with a small error bounded by the sum of ()65|1 and (j67j) . Using (j63|) . we finally find 

^«(n) ~ ne-''/^Yo{x{n, /t)) + e-''/'F(fi:) Jo(x(n, /t)) , (69) 



for n ^ Ik I, where we have defined 



F{k) = ^ - 7 - In K . 



(70) 



It is possible to express F{k,) in terms of the incomplete Gamma function, but we 
do not need this here. As a check, one may note that the Wronskian of the Bessel 
functions corresponding to (|56|) is W{jQ{x{n, K,)),YQ{x{n, k,))) = 2/tt, which after the 
appropriate coordinate transformations confirms W f^{n) , ^^{n)) = 1. 
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6.5 Resolvent and density of states 

We can now compute the resolvent of A and the spectral measure (density of states) by 
following the usual Sturm-Liouville procedure. We introduce, for k G C+, lm(K) 7^ 0, 
the functions /^(n) and Qnin), satisfying A/^ = nf^, and Ag,^ = ng^^, and characterized 
by the property that satisfies the boundary condition (jlSj) . /k(1) = /k(0), while (7^ 
is normalizable at infinity. We can express and g^ as linear combinations of $k and 
\E'k and we normalize them such that the Wronskian is equal to one, 

W{g^J^) = l. (71) 

Using /k and the resolvent of A can be written as 

f^MqAfn) for 72 < m 
R,{n,m) = {A- ^i)-\n,m) = r ^ ' ^ ' " (72) 

\9K{n)f^{m) forn > m 
Using and (jlH|) . it is easy to see that is given by 

= + . (73) 
To determine g^-, we recall that the Bessel functions Jq and Yq behave at infinity as 

Jo(x) ~ ./Acos(x-^) (74) 

" " y Tlx ^ 4 



X— +00 



yo(a;) ~ J— sin(x-^). (75) 

' x^oc \ TXX ^ A' ^ ' 

Using (jUUj) and (jU^ . and a; = I^JnR^ we see that the normalizable solution, behaving 
as Jq[x) ± iY{)[x) for Im(/t) ^ 0, is given by 

g^ = C± [ire-^^^ ± T ^e-'^F(/^)$,] , (76) 

where C± is determined by the requirement on the Wronskian, 

1 = W{g,, /«) = C± [Kire^^ T tKe-^F{K) t i\ ■ (77) 

The density of states or spectral measure p{k) of A is determined from the jump 
of the imaginary part of the resolvent across the real n axis. In the case at hand, one 
finds in the usual way, 

lim ^{g^+^e - gn-ie) = P(k)/k , (78) 

which upon using (fTSj) and (fTBj) gives, for k > 0, 

^^'^^ ^ (K7re-«)2 + (l + /te-'^F(K))2 ' ^^^^ 
For illustrative purposes, we show a plot of p{k) in fig. IHl 
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p(/t) 




K 



Figure 6: The density of states (f?^ of the operator A flMj) . 



7 Estimating the decay rate 



Equipped with the analysis of the hnear operator A, we now estimate the decay rate 
in the approximation explained in section |31 



7.1 A harmonic source 

We first study the radiation problem (j33|) with a harmonic source. 

Xn + o;o'A„ + ^(AA)„ = '^5n,i ae'"* . (80) 

There are in principle two ways to treat this problem. One is to make an ansatz of 
the form A„(t) = e"^*/i„ and to solve for imposing outgoing boundary conditions 
at 72 — > oo. The other is to use the spectral decomposition of A„(t) in eigenmodes of 
A. As a check, we will compute with both methods. 

7.1.1 Diagonalization of A 
We introduce 

oo 

AK(t) = 5Z/.(^)A.(t), (81) 

n=l 
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where is given by ()73|). The inverse transformation involves the density of states, 

An(t) = J p{^)dK Un) Kit) . (82) 
Using eq. (j8T| in eq. (j80|) yields 

K+Q^ + u;iy^ = be'^\ (83) 
where b = ^aff^{l). This is a forced harmonic oscillator and can be solved in the usual 



way. Using (j82j) to plug the solution back into the expression for the energy, 

oo 

m = E[2^" + + ^^^'^^^ ' ^^^^ 

n=l 

one finds that at late times, E is proportional to t and the radiated power is given by 

P = hm ^ = ^ua'piK) (85) 

where k = |(c<j^ — uJq) has to be positive, and p^n) is given by (fTUj). The derivation of 
fj85|) is similar to that of Fermi's Golden Rule. 

7.1.2 Modification of boundary conditions 

If we put A„(t) = e^'^^fin into eq. (jSUI) . we obtain 

-CuVn + ^(^At)n + tUo/^n = ^(^n.lO . (86) 

Thus, for n 1, fi = {fi„) must be a linear combination of $k and \E'k, with k = 
^{uj^ — UJq). We require that /i„ be outgoing at infinity, which from the behavior of Jq 
and Yq means that /i„ must behave as Jo{x) + iYo{x). Hence, 

/i„ = D [7re-'^$,(n) + z^,(n) - ze-'^F(/t)<l>«(n)] , (87) 

where D is determined by the equation at n = 1, 

2 2 2 

-- 2(/i2 - /ii) + (^0 - = (/^i - /^o) = ^ a , (88) 

yielding 

a = D [tckc-^ - tKe-T^K) - z] . (89) 
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To finish the calculation of the energy output, one has to plug the solution into the 
expression for the energy currenf ,^ 

Toi = -^(n + l)(A„+i-A„)A„. (90) 

Using the asymptotics of and \1'k and of Jq and Yq, one obtains for the average 
radiated power 

P = Wi{n ^ oo) = ^^e-'^l^r , (91) 
which using (jH^ . reproduces exactly (|H3j). 

with ft = |(u;2 - ul) > 0, and F{k) is given by (f7n|l . 
7.2 The decay rate of the NC soliton 

As explained in section the total decay rate of the NC soliton is obtained by de- 
composing the oscillation Aq into Fourier modes, using (jU^ for each mode, and 
summing over Fourier modes, i.e., 

oo 

Ptot(«,^)= Yl ^( 

ma,am,0), (93) 

m=— oo 

where 

' •- —imat 



am = 7r e-"""*Ao(t) (94) 

are the Fourier coefficients of Aq. The precise values of the of course depend on 
the details of the potential and the exact value of the basic frequency a (which, if we 
enforce energy conservation, is not fixed), so that we need an estimate to make further 
progress on (|^. 

There are two regimes in which one can estimate the Fourier coefficients. The 
first is for almost harmonic oscillations, for which one can use perturbation theory 
to determine the am- This will be the description when, for late times, the zeroth 
oscillator is in a small vicinity of the stationary point, and the soliton has almost 
completely decayed. The second regime is defined hy a <^ uoq. It describes the initial 



We define Tqi by the requirement that it satisfy Too(n) = — (Toi(n) — Toi(n — 1)) on a solution of 



the equations of motion, where Too{n) = ^A^ + j{^n — A„_i)^ + ^oJ^Xn is the "energy density" 
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stages of the decay of the NC sohton, and it will be a consistent approximation if the 
energy loss per period is much smaller that the total energy in the zeroth oscillator, 
i.e., we require 

Ptot{a, 9) < aHo < cooHo . (95) 

In this regime, the coefficients ()94|) can be estimated as follows |52j . 

If the energy Hq is exactly equal to the critical value V{Xmax), there is not only the 
static solution, but also the separatrix solution that is not periodic, but rather takes 
an infinite amount of time to fall off and come back to the top of the potential. We 
call this separatrix X{t). 

For energies close to the critical point of V, the oscillation, which has a frequency 
a <^ ujo, is well approximated by the convolution of the separatrix with a sum of delta 
functions separated by 27i/a, which is a much longer time scale than the scale typical 
of A, i.e., 

oo 

Ao(t) ^ Kt- 27cm/a) . (96) 

m=— oo 

The Fourier transform of Aq is then given by the product of the Fourier transform of A 
with delta functions at ma, m G Z. In other words, 

a f°° - ^ 

ttm^— dt e-"""*A(t) = a\{ma) (97) 
27r y_oo 

where A is the Fourier transform of the separatrix. 

Plugging ^7\\ in (jHSl), and using that p(k) is finite everywhere to replace the sum 
over Fourier modes by an integral, we obtain the final expression for the decay rate of 
the NC soliton, 



TT 



oo 

2 , 



Ptot{a,e)^—a dujuj\X{uj)\'p{K), (98) 

in which k = ^{uj"^ — uJq). 

For some special potentials, one can explicitly compute the separatrix and its 
Fourier transform. For instance, for the quartic potential that we used in the numerics, 
y($) = 4$2($ - 2)2, one has^ 



Another simple example that can be solved explicitly is the cubic potential. 
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and 



A(a) 



4^2 cosh 



8 



(100) 



For illustration, we show in fig. [7| the numerical evaluation of fj98p for this separatrix. 
We see that the value of 9 at which the energy loss per period becomes appreciable is 
around 6 ~ 0.2, in rough agreement with the estimate (jSHl)- 



P,,,{d,a)/aV{\ 
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Figure 7: The energy loss per period of Aq at initial stages of the NC soliton decay, as 
a function of 9. The potential is of the quartic form (pT|l. 



8 Conclusions 

In this paper, we have studied some properties of the dynamical system ()11|). which 
is (the radial part of) the classical equations of motions for a scalar noncommutative 
field in 2 + 1 dimensions. We have in particular analyzed the long time evolution of 
initial conditions of the form (II) corresponding to unstable noncommutative solitons. 
We have obtained numerical and analytical evidence that the energy is dissipated 
away to infinity in the form of classical radiation. For large, but finite values of the 
noncommutativity parameter 6', the decay rate of the soliton is exponentially small in 
9. Our main technical result is the explicit diagonalization of the linear problem. 
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Our paper will end here, but we would like to mention a number of open problems 
that are worthwhile of further investigation, and to which we hope to return in the 
near future. 

In the context of NC field theory, there are several other unstable solitons whose 
decays one might want to study. We have already mentioned the scalar soliton decays 
of type (12) and (13). Decays of the type (12) are particularly interesting because they 
have to proceed through non rotationally invariant configurations [37|. Solution of 
this problem will involve solving the linear problem for nonzero angular momentum 
{i.e., finding the "discrete Bessel functions" for degree / ^ 0, the function and 
corresponding to / = 0). Furthermore, it is not a priori clear into what decay products 
such a soliton will decay. For example, the soliton corresponding to Aq'^^ = = 
■ ■ ■ = Am/at, A^ ^ = Xmin could decay into either the stable soliton corresponding to 
Ao°^ = Xmin, = ^MiN > 0) or iuto the global vacuum. An additional question 
is whether there exist rotationally asymmetric stable solitons, and if so, whether they 
can be formed as decay products of the above process. It is natural to ask whether 
there is a simple criterion to decide this question. 

NC gauge theories, with or without couplings to scalar fields, also admit solitons, 
which in fact have a more direct interpretation in terms of D-branes in string theory 
(see for a review). Some of these solitons are unstable and hence can decay. Our 
solution of the linear problem is directly relevant in this context, but the nonlinearities 
are rather different, leading to a markedly different qualitative behavior j35j. 

Also from the point of view of string theory, it would be interesting to compare 
our results with the qualitative properties of the decay of the lump in 0^ theory in 
1 + 1 dimensions, which is another natural toy model for tachyon condensation. The 
fluctuations around the lump and the expression of the stable vacuum in terms of the 
fluctuations were studied in [^HJ. However, it seems to be unclear whether in the time 
evolution, the field eventually settles down in the stable vacuum, and if so, how the 
energy is dissipated pPj. 

Probably the most immediate problem from a mathematical point of view is to 
substantiate our qualitative discussion with rigorous estimates of the decay. As we 
have mentioned, our system is similar to systems exhibiting the discrete breather phe- 
nomenon One can also expect a connection to the general theory of resonances 
in nonlinear field equations and their radiative decay, studied, for example, in [H31 l54j . 
This is the general framework in which we expect a more rigorous treatment of our 
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problem. 
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